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The phase diagram of unzipping of an adsorbed directed polymer in two dimensions in a random 
medium has been determined. Both the hard wall and the soft wall cases are considered. Exact 
solutions for the pure problem with different affinities on the two sides are given. The results 
obtained by the numerical procedure adopted here are shown to agree with the exact results for the 
pure case. The characteristic exponents for unzipping for the random problem are different from 
the pure case. The distribution functions for the unzipped length, first bubble and the spacer are 
determined. 
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I. INTRODUCTION 

The study of polymer adsorption on a surface is impor- 
tant because of many applications such as in lubrication, 
adhesion, surface protection, coating of surfaces, wetting, 
vortex lines, and biology. Of particular interest is the 
adsorption-desorption transition of a polymer, similar to 
the denaturation transition of a double stranded DNA 
molecule. The adsorbed phase, dominated by the attrac- 
tion with the wall, wins at low temperatures whereas at 
high temperatures the desorbed phase prevails. Although 
the temperature driven adsorption-desorption transition 
dates back to sixties Q, 0, 0, 0, 13 , the advent of micro- 
manipulation techniques such as optical tweezers, AFM, 
etc. 0, regenerated interest in it because of the pos- 
sibility of exploring single molecules especially under an 
external force. 

The existence of a force induced unzipping transi- 
tion was established theoretically for DNA type dou- 
ble stranded molecules in Ref. pj- Many aspects 
of the unzip ping tr ansition have now been elucidated 
I^ IIclll^fiallaTT^ I. including an experimental phase di- 
agram |15| . The basic feature of the unzipping transition 
is that, at a fixed temperature T, the polymers stay in the 
bound or localized phase if the magnitude of the pulling 
force, g, is less than the critical force g c (T), but they can 
be unbound or unzipped if the pulling force exceeds the 
critical force. Moreover, the transition is described by 
a set of characteristic exponents. Since this transition 
involves a competition between a binding potential and 
the stretching by the force, such an unzipping transition 
is also possible for an adsorbed polymer Il3 ■ 

Apart from thermal fluctuations and force, disorder 
(or impurities) in the medium may also induce the 
adsorption-desorption transition even below the desorp- 
tion transition temperature. This phenomenon is known 
as a depinning transition. The depinning transition for a 
directed polymer has been studied in the past because of 
its importance in understanding fluctuations in domain 
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walls in random magnet, pinning of magnetic flux lines 
in dirty superconductors having columnar defects and in 
particular as a simpler siblin g of the more complex poly- 
mer adsorption problem [Lsl Il9l l2p| | . In a homogeneous 
medium, in absence of any attraction with the wall, a 
polymer behaves like a free chain controlled by entropy. 
But with quenched impurities, there may be one or more 
ground state configurations |2l|. As a result, a polymer 
may swell over a free chain to take advantage of the at- 
tractive locations in the medium. This swelling effect 
is pronounced in lower dimensions and would affect the 
response of the polymer to any applied force. More im- 
portantly, because the impurities are quenched (no ther- 
malization allowed), the behavior of a polymer is sample 
dependent. An averaging over the samples often yields 
well-behaved thermodynamic quantities but still sample 
variations may reveal drastic differences from the average 
behavior. Consequently the unzipping phenomenon in a 
random medium ( "dirty environment" ) is expected to be 
different from that in a homogeneous medium because, 
in the random case, it is not just pulling off from the wall 
but also pulling away from local pockets of favorable con- 
figurations. 

In this paper our aim is to study the phase diagram 
of a directed polymer which is adsorbed on a wall in a 
random environment and pulled by a force at one end to 
unzip it off the wall, keeping the other end fixed. Though 
the unzipping transition in a pure system is first order, 
randomness may change the order of the transition. This 
change will be reflected in the exponents that describe the 
transition. Therefore, in addition to the overall phase 
diagram we would also like to determine the exponents 
in this paper. 

The disorder in the problem can be introduced in sev- 
eral ways as shown in Fig. ^ These possibilities are as 
follows, (i) The polymer sequence is homogeneous, and 
the randomness is introduced on the lattice site such that 
the ground state of the polymer (on the wall) remains 
unaffected as if it is in pure environment. In this case 
the binding from the wall is attractive and constant, (ii) 
The randomness is introduced in the polymer sequence 
but the environment is kept pure. The binding from some 
portion of the wall, depending on the strength of random- 
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FIG. 1: Schematic diagram of a directed polymer in a random 
medium with an attractive wall. By construction, one direc- 
tion along the polymer is special. The wandering is in the 
transverse d directions and therefore it is d + 1 dimensional 
model. Here d — 1. There is a line at x — (a wall or an inter- 
face) that attracts the polymer, (a) Hard wall: polymer is not 
allowed in the region x < 0. (b) Soft wall: Polymer is allowed 
in the whole region. There may be an extra potential V(> 0) 
on one side, say x < 0. (c) Homo-polymer in random medium. 
There is a unique (sample independent) ground state for this 
case, (d) Hetro-polymer in pure medium. The ground state 
is sample dependent. In (c) and (d), randomness is in the 
shaded regions. 

ness, can be attractive or repulsive. The ground state in 
this case is sample dependent (iii) The randomness is 
put at each lattice point including the wall. For strong 
enough disorder some sites of the wall may become re- 
pulsive. In this case also the ground state is sample de- 
pendent. Of these, item (ii) of heterogeneous sequence 
is of importance because of its similarity with the DNA 
unzipping problem. In this paper we restrict ourselves to 
case (iii) only, especially because of the known behavior 
of a polymer in a random or dirty medium. 

The thermodynamic or bulk behavior is described by 
temperature and one of the two conjugate variables force 
or distance. This results in two different ensembles, 
namely the fixed force and the fixed distance ensemble. 
For the pure case the ensembles are equivalent when the 
polymer is pulled at the end. But this is not always the 
case. It is shown recently by Kapri et. al 0] in the 
context of DNA unzipping that the ensembles arc not 
equivalent when the pulling force is applied at any frac- 
tion s (0 < s < 1) on the double stranded DNA. In this 
respect pulling at the end is a special case in which re- 
sults happened to be ensemble independent. We restrict 
ourselves to the force at the end case only and again 
show that the pure phase diagram can be obtained ex- 
actly from both the ensembles. However there are subtle 
differences and such differences play a significant role in 
the random case. 

Apart from exact solutions for the pure case, our ap- 
proach involves the use of exact transfer matrix for a 



polymer on a two dimensional square lattice. The choice 
of two dimensions is partly influenced by the known re- 
sults of strong effect of disorder in lower dimensions. The 
paper is organized as follows : The model is described 
in section II. The quantities of interest are also defined 
here. In section III we discuss the adsorption-desorption 
transition in absence of randomness. The exact phase 
diagram and a comparison with exact transfer matrix re- 
sults can be found here. We also explain the use of finite 
size scaling to extract the critical value and the expo- 
nents from the results of extensibility for finite chains. 
A general case of the adsorbing wall separating media 
of two different types is also discussed. Section IV is 
devoted to the case of a polymer in a random environ- 
ment. We determine the phase diagram and the changes 
in the exponents due to randomness. These are based 
on the results of the transfer matrix for finite chains and 
are therefore exact for each sample. The statistical error 
comes only because of averaging over a finite sample size. 
For the data collapse in finite size scaling analysis, the 
Bhattacharjee-Seno procedure has been adopted |22j . 



II. MODEL 

We model a polymer by a directed random-walk in a 
two dimensional [d = 1 + 1) square lattice directed along 
the diagonal of the square. There is an attractive wall 
on the diagonal of the square (a; = 0) which tends to 
suppress wandering because the polymer gains an energy 
— e (e > 0) each time it is on the wall [23]. Depending 
on the ensemble one is working with, either a fixed force, 
g, is applied at the end monomer (fixed force ensemble) 
or the distance, x, of the end monomer from the wall is 
kept fixed (fixed distance ensemble). The other end is 
always kept anchored. Two types of walls are possible, 
namely (i) a hard wall if the polymer can not cross the 
wall, and (ii) a soft wall if the polymer can cross the 
wall. For the random case, there is an additional random 
energy at each lattice point. In a continuum notation, the 
Hamiltonian may be written as 

^H^dD + jf * dzV - {r{z}} 

f N f N dr 

+ dz rj(r(z),z) -g- / — dz (I) 
Jo Jo dz 

where r(z) is the d-dimensional co-ordinate at monomer 
index z, V w (r) is the potential due to the wall, r\ is the 
random energy and g is the applied force at the end point 
z = N . In this paper we consider a discrete version (Fig. 
^) of this model. 

A more general case can be considered where the poly- 
mer experiences different potentials on the two sides of 
the wall, as, e.g., may happen for adsorption at the in- 
terface of two immiscible liquids with different affinities 
for the polymers. A hard wall case may then correspond 
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to a solid-liquid interface. The general case is considered 
in the next section and in detail in Appendix A. 



A. Pure case 

The partition function for the pure case is obtained by 
a recursion relation 



?iv+i(x) = ^2 Z N (x+j) [l + (e 
j=±i 



0e 



l)k,o] (2) 



where Zm{x) is the canonical partition function of the 
polymer of length N whose 7V th monomer is at a distance 
x from the wall, the zeroth monomer being at x = 0, and 
(3 is the inverse temperature l/T in units of the Boltz- 
mann constant fee = 1. The temperature dependence 
of the partition function has not been shown explicitly. 
The initial condition for the above recursion relation of 
the partition function is 



Z (x) = eP*8 t 



x,0- 



(3) 



The non-crossing constraint of the polymer at the hard- 
wall is taken care of by reassigning 



Zj(x)=0, Vx < 0, 



(hardwall) 



(4) 



after each step j. 

The canonical partition function with a force, g, acting 
at one end keeping the other end fixed is then calculated 
by summing over all the allowed configurations of the 
walks on the lattice. 



z N{g) = ^2 Z N (x)e 



(3gx 



(5) 



where e@ gx is the Boltzmann weight due to force g. The 
thermodynamic properties in a given ensemble are ob- 
tained from the free energy, 



F N (T,x) 



-Tin Zn(x), 
-TlnZ N (g), 



(fixed distance) (6a) 
(fixed force). (6b) 



is applied is the quantity of interest and can be obtained 
by 



Ex Z N (x)ePB* ' 



(8) 



The finite size effects in fixed distance ensemble are 
weak which gives the privilege of working with chain of 
shorter length and still getting the exact phase boundary. 
In contrast, the fixed force ensemble exhibits strong finite 
size effects and finite size scaling has to be used to get 
the phase boundary. 

The response of the polymer is characterized by a sus- 
ceptibility like quantity to be called isothermal extensi- 
bility defined as 



XT 



d(x) 
dg 



(9) 



This extensibility can be expressed in terms of the posi- 
tion fluctuation in the fixed force ensemble as 



1 



XT = 



k B T 



Ix 2 ) - (x) 2 



(10) 



This relation is useful for numerical computation. We 
determine the isotherms g vs (x) or (g) vs x (depending 
on the ensemble). 



C. Polymer in a random environment 

In the random environment the recursion relation has 
also the contribution due to the randomness. The recur- 
sion relations in the fixed distance and the fixed force 
ensemble respectively read as 



J N+ 



l(*) 



]T Z { N a} (x + j) [l + (e*- 1)^,0] 
j=±i 

xe -^(x,JV+l) j (na) 



where t)(x, N) is the random energy at site (x, N), and 



ZP(g) = ^eP-zPiz) 



(lib) 



B. Ensembles 

For the fixed distance ensemble, the distance, x, of 
the last ( N th ) monomer from the wall is kept fixed. 
The quantity of interest, the average force required to 
maintain this distance is given by 



(fj) 



dF(T, x) 
dx 



(k B = 1) 



(7) 



and is obtained by the finite difference in free energy F(x) 
of Eq. US] for the lattice problem. 

In the fixed force ensemble, the average distance of the 
last monomer from the wall, (x), where a fixed force, g 7 



The superscript a in the above equations denotes a par- 
ticular realization. For every sample, Eq. [3] and 0] remain 
valid so that these recursion relations can be used both 
for soft- wall and hard- wall cases. 

The random energies are independently distributed 
with no spatial correlation (white noise), drawn from a 
uniform deviate of zero mean (^(x, r)) = 0, width A, and 



[n(x, TMa/yjU = — 6(x - x') 5{r - r'). (12) 

Throughout this paper [...]di s denotes the average over 
various realizations and (...) denotes the thermal aver- 
aging. 
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The thermodynamic quantities with quenched ran- 
domness come from the average free energy of the ap- 
propriate ensemble, i.e. the average of F^ a} (x) for the 

fixed distance ensemble, and pffl (g) for the fixed force 
ensemble. 

With e = 0, the randomness leads to a swelling of the 
polymer chain as measured by the transverse size 



[< x 



N >]dis 



N 



2// 



(13) 




where v = 2/3 in two (1+1) dimensions [21|. For a pure 
system, v = 1/2. 



Bubble 



FIG. 2: Schematic diagram showing the first bubble of length 
mj , spacer of length m 3 and the unzipped segment of length 
m when a pulling force g is applied at one end of the polymer 
in the random medium. 

It is apparent that the critical force needed to unzip 
a polymer has to be sample dependent. Even in a fixed 
force ensemble, with g at its critical value, the whole 
chain may not be unzipped in all samples. Important 
quantities to study are then the unzipped length m mea- 
sured from the open end to the first contact ("Y-fork"), 
the length m s of the completely zipped region after that, 
to be called a spacer, and then the length m& of the first 
bubble. A bubble is defined as an unzipped region sep- 
arated by two contacts or bound monomers. These are 
defined in Fig. [21 Apart from the phase diagram ob- 
tained by studying the response functions averaged over 
samples, the distribution functions of the various lengths 
over samples are determined. 



limit. For any V ^ 0, there is a temperature driven 
desorption transition with T c — ► oo as V — > 0. The 
unzipping transition is given by (u = e~^ e and v = e~@ v ) 



g c {T) = T\n 



l + v(l-2u) 
u(l — uv) 



- 1 



(14a) 



or, equivalcntly, 



g c (T) = 2Tcosh~ 



1-1- 



2it(l — vu) 
1 + v(l - 2u) 



(14b) 

To be noted is the small region at intermediate temper- 
atures where a reentrance can be observed. The details 
are given in the Appendix A. 

The unzipping on the wrong side, i.e. with x < where 
V > 0, has a different phase boundary than Eqs. I14al and 
I14bl and is given by 



9c(T) 



-2Tcosh" 



1 



1 



1 



2ti(l - vu) 
l + v{l - 2u) 



(15) 

The phase diagram is not symmetric around g = 0, ex- 
cept for V = 0. The critical force need not tend to 
zero at the transition temperature to unzip on the wrong 
side. At T = 0, the critical force for this model is 
5c(0) = — 2(V + e/2). More details of the phase diagram 
for V ^ will be discussed elsewhere. 

Fig. |31 shows the phase diagram for general V, includ- 
ing the soft- wall and the hard- wall cases. The lines arc 
from exact results and the symbols from the numerical 
analysis presented below. We show in Figs. 01 and 
the final phase diagrams obtained by the same numeri- 
cal procedure for the random system. These constitute a 
part of the major findings of this paper. 



B. Units and dimensions 



III. EXACT SOLUTION: PURE CASE 
A. phase boundary 

The recursion relations, Eqs. [2Jand[I3 for the pure case, 
can be analyzed exactly. A more general case where the 
polymer is adsorbed at the interface of two different types 
of media (e.g. immiscible liquids with different affinity for 
the polymer) is considered |16| and a generating function 
based derivation is given in Appendix A. Also we show 
how the phase diagram can be obtained from a fixed dis- 
tance ensemble. 

The difference in affinity is modeled by a potential V > 
for x < 0. V < oo would correspond to two immiscible 
liquids while the hardwall V — > oo corresponds to an 
interface between a liquid and an impenetrable solid. The 
homogeneous soft wall case corresponding to the V — > 



Throughout the paper we choose dimensionless quan- 
tities : Boltzmann constant ks is chosen such that 
e/fe_B = 1, and in this unit, Ub = 1 with temperature 
T as dimensionless. Distances are measured in units of 
lattice spacing a — 1 (along the diagonal) . The quantities 
like force (g), width (A) and potential V have dimension 
of energy and are measured in units of e. All the plots 
shown in this paper use the dimensionless forms where 
g — > g/e, T is A^T/e, V — > V/e, A — ► A/e and x — ► x/a. 
In case there is no adsorbtion energy, the relevant pa- 
rameter is A/ksT and the choice of e does not matter. 

C. Isotherms 

For the unzipping transition, (x) vs g or (g) vs x 
isotherms (depending on the ensemble) are of interest. 
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FIG. 3: Dimensionless force (g) versus temperature (T) phase 
diagram for an adsorbed polymer. Each lattice site on one 
side (x < 0) of the wall has a repulsive potential of strength 
V > 0. Positive g corresponds to case of pulling the polymer 
away from the wall on favorable side while the negative g (note 
the different scale) is for pulling on the wrong side. The two 
limiting cases are V — » (softwall) and V — > oo (hardwall). 
The points are from numerics and the lines are the exact 
analytical results. Note the reentrance on the positive side 
for V — 0.67 (shown in inset). 




FIG. 5: Phase surface for (a) hardwall and (b) softwall in 
the presence of disorder. 
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FIG. 4: Critical force vs disorder strength for both the hard- 
wall and the softwall cases at f3e = 15, i.e. (T = 0.067). The 
thin slice of reentrance is shown in the inset. 



The phase diagram can be mapped out from the special 
features of these isotherms. The procedure is explained 
for the pure case and the numerical results are compared 
with the exact results. This procedure is to be adopted 
for the random case. 

Fig. shows the isotherm for the pure case. For finite 
N, these are obtained by iterating the recursion relation 



("transfer matrix" approach) to calculate the partition 
function and then the extension or the force depending 
on the ensemble used. Except for a region near x = 0, 
the two are identical (within numerical accuracy). The 
difference in the short distance region is expected because 
x in fixed distance ensemble has to be an integer but 
< x > in the fixed force ensemble need not be |2J| . 

The flat region in the g vs x isotherm signals a co- 
existence between a zipped and an unzipped (desorbed) 
phase. This is for a first order transition. The critical 
value g c (T) is obtained from the intersection of the flat 
part of the g — x isotherm with the isotherm of the pure 
phases. For a higher order transition, there may not be 
any flat region but some singular feature would persist. 
In either cases singularity on the isotherm gets rounded 
for finite N, and a finite size scaling analysis can be done 
to extract the critical value. The advantage of using this 
finite size scaling is that it does not rely on the order of 
the transition. This method remains valid even if there 
is no flat part of the isotherm (signaling a higher order 
transition). The exact solution suggests a scaling form 
for extensibility \, 



X = N d *X((g~g c )N*) 



(16) 



with two characteristic exponents d x and <j). As we argue 
in Appendix A, the scaling variable for the pure problem 
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FIG. 6: Force distance isotherms in (a) fixed force ensemble 
(b) fixed distance ensemble, for a polymer of length N = 2000 
at @e = 1.5 (T = 0.067). Each lattice site on one side of 
the wall contains a repulsive potential of strength "V " and 
polymer is pulled away from the wall on other side. The inset 
in (a) shows the knees which developed in fixed force ensemble 
for low x but absent in fixed distance ensemble. 



is N(g — g c ) so that 0=1. 

The scaling variable in Eq. 1161 can be taken to imply 
that the length m of the unzipped region (see Fig. [21 
depends on the force as a power law, 

m-lff-ffel- 1 ^. (17) 

This length m has been used in studying the dynamics 
of unzipping of DNA ^lj . 

If we demand extensivity for x for N — » oo, then 



X(x) 



-(d x -i)/0 



X/N 



This implies 



9 c | 



-{d x -i)/<f, 



(18) 



A similar scaling form can be written for the specific 
heat (per monomer). At a fixed force, the temperature 
dependence of the specific heat is 



c g (T) = N a ^y((T - T e (g))N+*), 



(19) 



which defines the exponent <pt in addition to the usual 
bulk specific heat exponent a. To recover the bulk behav- 
ior, y(x) ~| x \~ a for large x. For a first order transition, 
a = 1. 



The unzipping transition is characterized by the four 
exponents, d x ,<p,a,<pf These are known from the ex- 
act solution and we compare with the numerical results. 
These are also the quantities needed for the random sys- 
tem. 



IV. TRANSFER MATRIX: NUMERICAL 
RESULTS FOR THE PURE CASE 

We use the transfer matrix technique in both the en- 
sembles to calculate the exact partition function numer- 
ically from the recursion relations in Eqs. [21 and 

In Fig. [7| an example of data collapse for extensibility 
X for the hard-wall case is shown. This is at (3e = 15 
(T = 0.067) for N = 1000, 1500, and 2000. The critical 
exponents suggested by collapse are 



and 



1. 



(20) 



consistent with the exact solution in Appendix A. This 
implies 

x/N-lg-gc]- 1 , and m~\g-g c \- 1 . (21) 

In the thermodynamic limit, both the ensembles give 
the same phase boundary. By the data collapse method 
for the extensibility, the critical force g c (T) can be de- 
termined for a given T. Once the critical force is known, 
the size (N) dependence of the average distance (xj) of 
the i th monomer for forces below and above it can be 
determined. We have checked that for g below g c , (a;,) 
decreases very rapidly with N showing that the polymer 
is adsorbed on the wall whereas for g slightly above g c 
the polymer is desorbed from the wall with < Xi >cx i. 

Also from the same calculation the specific heat has 
been computed. A comparison of the specific heat is 
shown in Fig. [H] The finite size scaling variable that one 
gets from the specific heat is of the form (T - T c {g))N^ 
with <$> t = 1. Though a good collapse is obtained for 
a = 0.93, but the error bar does not preclude a = 1 
(shown in Fig. [SJa), especially if there are finite correc- 
tions. The development of the S- function peak (a = 1) 
at the transition point is a vindication of the first order 
nature of the unzipping transition. 

The critical values obtained from data collapse of the 
numerical results are used to construct the phase dia- 
gram. These are shown in Fig. The numerical results 
(points) are in good agreement with the analytical results 
(lines) . 



V. RANDOM CASE: PHASE DIAGRAM 

To get the phase diagram we need to calculate the crit- 
ical force, g c (T), at different disorder strengths and tem- 
perature keeping the binding to the wall constant. We 
use the finite size scaling of the correlation function for 
that purpose. For a disorder problem, depending on the 
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FIG. 7: Data collapse of extensibility at fie = 15 (T = 0.067) 
for the pure case (a) softwall (b) hardwall. The lengths of 
the chain used are N - 1000, 1500, 2000. For both cases, the 
exponents are d x — 2 and cj> = 1 



way disorder averaging is done, two different correlation 
functions can be defined. 



and 



Xdis = \{x 2 ) - (x) 2 ]dis 



Cdis = [(X 2 )]dis ~ [{X}]l h 



(22a) 



(22b) 



For Eq. I22al the disorder averaging is done at the re- 
sponse function level of an individual sample whereas for 
Eq. I22bl the disorder averaging is done at the moment 
level. The critical exponents for these two different cor- 
relation functions may not be the same. For the problem 
in hand, correlation functions scale as 



and 



■'dis 



(23a) 



(23b) 

where d x and dc are the critical exponents for Xdis and 
Cdis respectively with <fi describing the finite size behavior 
of the force. For the random problem, (f> is expected to 
be different from the pure case of Eq. ^| 



N dc W {{g - g c )N*) 
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FIG. 8: Specific heat up to a multiplicative constant for 
g — 1.0. (a) There is a peak developing as N increases, over 
and above a discontinuity at T c (g) = 0.87204. (b, c) Data 
collapse at T c (g) for chains of length N = 1500,2000,2500, 
with cj>t = 1 and acj)t = 0.93 in (b) but acj>t = 1 in (c). 



In the fixed distance ensemble the force, g^(x, T), re- 
quired to maintain a distance x from the wall is computed 
for each realization by taking the finite differences in free 
energy F^{x,N) = -k B TlogZ^{x,N), which on av- 
eraging over samples give the average force, (g(x,T)), 
required to maintain the distance a;. 

(«?(*, T)) = [F { N a} (T,x + a)-F { N a} (T,x)} dis /a (24) 

where a is the lattice constant (in our case a = 1 ). 
This has also been used to generate the phase diagram 
or crosscheck the results. 

The data collapse for both the correlation functions are 
shown in fig. The exponents for the hard-wall case, as 
suggested by the collapse, are 



= 0.5, d x = 1.5, and d c = 2.0. (25) 

For both the correlation functions we obtain the same 
g c value suggesting that any of them can be used to cal- 
culate the critical force. It is found that the correlation 
function at the moment level (Eq. I23b|) behaves more 
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FIG. 9: Data collapse for two different correlation functions 
for the random system, (a), (b) are for hardwall and (c), 
(d) are for softwall. The chain length N = 1000, 1500 and 
2000. In (a) and (c) the disorder averaging is done at the 
response function level of an individual sample (by definition 
it is the extensibility upto a factor of T _1 ). The exponents 
are d x = 1.5, 4> = 0.5. In (b) and (d) averaging is done at 
moment level. The exponents are dc = 2.0, (f> = 0.5. For all 
these correlation functions, the disorder strength A = 2/3 and 
averaging is done over 10 5 realizations at j3e = 15 (T = 0.067). 



smoothly than the response function (Eq. I23a|) and so 
Eq. 123 bl can be used to calculate g c , with same accuracy 
for less number of realizations. 

The phase diagram obtained by this determination 
of g c (T) for various T and the dimensionless disorder 
strength A is shown in Figs. 0]and[S| For small enough 
A and also small T, the random energies at the off- wall 
sites are not energetically favourable to delocalize the 
polymer. But still there are pockets or sites with lo- 
cally favourable energies. The pulling force then has to 
pull the polymer out of the wall and the local pockets 
signalling an increase in the critical force until entropic 
effects start dominating or the disorder favors a delocal- 
ized state. Such a thin slice of re-entrance for small A is 
shown in the inset in Fig. 21 The full three dimensional 
surface in the g - T - A space is shown in Fig. [SJfor the 
hard and the soft wall cases separately. 

From the determined values of the exponents, we find 
that the quenched averaged extensibility x, which has to 
be extensive, should have the same behavior as the pure 
case, namely, 



X/N 



and 



\9-9c\ 



(26) 



It is to be noted that for heterogeneous chain this ex- 
tensibility has been predicted to be \g — g c \~ 2 , though 
the finite size exponent </> is 1/2. In other words, the 
average behavior of the extensibility remains pure-like 
though the size dependence or the scaling variable is like 
a heterogeneous chain case. This is one of the important 
conclusions of this paper. For a soft wall shown 
in Fig. the exponents are the same as in the hard wall 
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FIG. 10: Specific heat for the random system at (a) g = (b) 
g = 1. There is a finite size dependence which however does 
not scale to yield a good data collapse. The averaging is done 
over 10 4 samples. 



case. 



The zero force specific heat is shown in Fig. llUfa . Un- 
like the pure case, the discontinuity at the transition tem- 
perature is not clearly visible. It is interesting to note 
that all the specific heat curves for different values of N 
seem to cross at the same temperature. Such crossings 
generally imply existence of a certain "universal" scale 
|25j . In the present case this scale could be due to the 
proximity of the phase transition point but its identifica- 
tion remains to be done. There is a size dependence of 
the specific heat near the peak as N increases. This may 
be a signal of a very weak divergence of specific heat or 
even the formation of a cusp (a < 0). Fig. HHId shows the 
specific heat vs T for g = 1.0. Here also we see a growth 
in the peak, but it is not as sharp as in the pure case 
(Fig. EJ. No good collapse can be obtained. Whether a 
delta peak is forming as N — > oo is not obvious. Longer 
chains need to be studied to sort out these issues. 
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VI. SAMPLE FLUCTUATIONS 
A. Response 

To study the response of the force, the partition func- 
tion is computed for polymer of lengths up to N = 2048 
at a fixed temperature T, binding to the wall e and 
disorder strength A. The average distance of the last 
monomer, X, where a fixed force g, is applied (fixed force 
ensemble) is 



where 



P^{x) = 



J N 



(x) 



E*4 a} (*) 



(27) 



(28) 



is the probability of the end monomer of sample a being 
at a distance x from the wall. 
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FIG. 11: (a) Re-scaled force gN 1 ' 3 versus re-scaled average 
distance (x}/N 2 ^ 3 , from the wall for four different samples of 
size N = 1024 and N = 2048 at pe = 15 (T = 0.067) for the 
disorder strength A = 2 when there is binding to the wall (b) 
Extensibility for the same samples. 

Fig. ^2 shows the re-scaled force gN 1 / 3 , versus re- 
scalcd separation, (x)/N 2 / 3 from the wall fFig lllf a'O and 



the extensibility fFig lllf bll for four different samples of 
lengths N = 1024 and N = 2048 at (3e = 15 (T = 0.067) 
for disorder strength A = 2 when there is a soft- wall. 
The large plateau for small force shows that the energy 
gain from the wall favors the adsorption of polymer and 
some critical force, which depends on sample, is needed 
to desorb it. Above this critical force the polymer does 
not see the wall but because of the randomness it can still 
get trapped to some attractive sites. The distance from 
the wall does not increase when the force is increased 
by small amount. The response would then just be the 
thermal width of the probability distribution which at 
low temperatures is very narrow and independent of N. 
As a result small plateaus develop in the isotherm for 
any specific sample as shown in the plot 0, |2^| . But be- 
cause of the possibility of nearly degenerate but spatially 
separated states, a fixed force of right amount may pro- 
duce a large displacement. Consequently, sharp jumps 
appear between plateaus in the isotherms. The corre- 
sponding extensibility of a sample shows spikes (delta 
function peaks) where there are jumps in the force ver- 
sus distance curves and it is zero within plateaus. Note 
the ensemble difference here. Whereas a small change in 
g (in the fixed force ensemble) can produce a large change 
in x, in the fixed distance ensemble, a small change in x 
would just explore the neighborhood of the preferred con- 
figuration. The spikes of Fig. II If remain undetected. 
On averaging over samples, we need to get back exten- 
sivity for the extensibility \- This is possible because of 
rare samples with probability ~ N~ 2lJ+l having degen- 
erate but spatially separated states (allowing unzipping 
at zero- force) plj . 

The plot between re-scalcd force and re-scaled distance 
from the wall and the corresponding extensibility but 
with e = 0, for the same parameters is shown in Fig. 
IT21 The average response of the force is linear at small 
forces 

It is interesting to calculate the probability distribution 
for average height of steps, Ax. Fig. ^1 shows a log - log 
plot of P(Ax) versus the step height Ax. Here we have 
assumed that the steps are independently distributed. 
These are taken over 10 5 samples at [3e = 15 (T = 0.067) 
and disorder strength A — 2. The straight line is the 
best fit which suggests the probability distribution of the 
form 



P(Ax) ~ (Ax)' 



(29) 



with p = 1. 

In Fig. ^]we have plotted again the re-scaled separa- 
tion with re-scaled force for both binding and no binding 
to the wall but by keeping the underlying disorder same 
at a disorder strength (A = 2) and f3t = 15 (T = 0.067) 
for N = 1024. The steps are for a single sample and 
the smooth curves are for averages over 10 4 samples. For 
large values of force both the curves merge showing that 
the effect of the binding is negligible at large force. 



10 




FIG. 12: (a) Re-scaled force gN 1 ^ 3 versus re-scaled average 
distance (x)/N 2 ^ 3 , from the wall for four different samples of 
size N = 1024 and N = 2048 at /3e = 15 (T = 0.067) for the 
disorder strength A = 2 and there is no binding to the wall 
(e = 0). (b) Extensibility for the same samples. 
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FIG. 13: Log log plot of the probablity distribution P({Ax)) 
of the rescaled average step height (Ax). This is for a polymer 
of length N = 1024 in a fixed force ensemble at f3e — 15 (T = 
0.067), disorder strength A = 2 and there is no binding to 
the wall (e = 0). 10° samples used to generate the statistics. 
The slope of the fitted line is 1. 
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FIG. 14: (a) Re-scaled force gN 1 ^ 3 versus re-scaled average 
distance (x)/N 2 ' 3 , from the wall for a polymer of length N = 
2048 at (3e = 15 (T = 0.067). The disorder is kept same (at 
a disorder strength A = 2) for both binding and no binding 
(e = 0) to the wall. The steps are for a single sample and 
smooth curves are after averaging over 10 4 samples. 
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FIG. 15: Data collapse for the probability distribution of 
pseudo critical force P(go)/N 3 versus (g c - g )N<l>* for the 
hardwall case for N = 512, 1024. 
The averaging is done over 10 samples. 



d g = 0.41 and (f> g = 0.31. 



It is clear from Fig. 1111 that in the random environ- 
ment, the critical force needed to unzip a polymer de- 
pends strongly on sample. E.g., the displacement for a 
particular sample would depend on the deviation g — go, 
where go is the critical force for that sample. A random 
system may have another scale that shows the closeness 
of go to the average critical force g c . The scaling be- 
haviour of this scale can be obtained from the probability 
distribution of g or the finite size scaling of the distri- 
bution, P(g ) = N d ^(\g - ffcl^ 03 ) that defines d g and 
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We take the first jump from the wide plateau around 
g = as the measure of the (pseudo) critical force, g ■ 
The scaled probability distribution obtained from the fre- 
quencies of go over 10 5 samples for N = 512 and 1024, 
is shown in Fig. The plot suggests d g = 0.41 and 
4> g = 0.31, which also implies a sharply peaked proba- 
bility distribution. The fact that <f> g ^ <f> (see Fig. |§J) is 
suggestive of different scaling of the two different scales. 



C. Bubble, Spacer and Unzipped segments 

In the random medium, the disorder controls the con- 
figuration of the polymer. As already said, the criti- 
cal force needed to unzip the polymer is sample depen- 
dent. It is interesting to study the distribution of various 
quantities such as bubble length, length of adsorbed or 
zipped and desorbed or unzipped segments with the ap- 
plied force. 

25 I 1 . . 1 1 1 

20 - 




FIG. 16: Average separation from the wall, (xt), of i th 
monomer of chain of length TV = 300 at f3e = 3 (T = 1/3) and 
disorder strength A = 5/3 when a stretching force g = 0.5 is 
applied at one end. Four different chain configurations are 
shown by the dotted and dashed lines. The solid line is when 
averaging is done over 90000 samples. 

Fig. ^j] shows the average separation from the hard- 
wall for each monomer of chain of length N = 300 at 
a temperature f3e = 3 (T = 1/3) and disorder strength 
A = 5/3 when a stretching force g = 0.5 is applied at 
the end. Four different chain configurations are shown 
by the dotted and dashed lines whereas the curve with 
solid line is the result of averaging over 90000 such chain 
configurations. Although the applied force here is below 
the critical force needed to desorb the polymer yet there 
are some samples which are in the desorbed phase. 

In this paper we concentrate only on the probability 
distribution of (i) the length of the desorbed or unzipped 
segment, (ii) the length of the adsorbed or zipped seg- 
ment between unzipped segment and the bubble (spacer) 
and (iii) the length of the first bubble. These quantities 
arc defined schematically in Fig. [2J 
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FIG. 17: Probability distribution of the unzipped length m. 
(a) Plot of P(m) vs m p , at /3e = 3 (T = 0.333) and disorder 
strength A = 5/3, when a stretching force g = 0.5 is applied 
at the end. The points are for N = 200 and 300 and the solid 
line is the fit using the functional form f(m) = Ae~ am . Here 
p = 0.51±0.01. (b) Plot of P(m) vs m at j3e = 3 (T = 1/3), a 
force g — g c = 0.82 is applied. This is the force on the phase 
boundary. The solid curve is the fit using the functional form 
f(m) = co/m b " + ci/(N — m) bl .To obtain the fit, we have 
excluded m = and m = N. For this plot c = 0.026 ±0.001, 
bo = 0.37 ± 0.0, ci = 0.028 ± 0.001 and bi = 0.71 ± 0.02 are 
the fitting parameters. N = 300 for both the figures. 



The probability distribution (over samples) for the un- 
zipped length P(m), for the length m of the chain which 
gets desorbed or unzipped at the phase boundary is 
shown in Fig. 1171 For force below the phase boundary, 
the probability distribution decays rapidly with length 
(Fig. IT7h ) but not so on the phase boundary (Fig. H7b). 
For g < g c {T), the probability distribution seems to be 
stretched-exponential, P(m) ~ exp(— a^Jm), but the dis- 
tribution for the phase boundary seems to admit power 
laws, especially near the peaks at the two extremes. Such 
a two peak structure is expected at a phase coexistence. 

Fig. ^| shows the probability distributions for the 
length of spacer (m s ) and the length of the first bub- 
ble (mt). Both the distributions at phase boundary are 
strctchcd-exponcntial P(m s ) ~ exp(— am s Cs ) with c s = 
0.72 ± 0.01 for spacer length and P(mi,) ~ cxp(— amb Cb ) 
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FIG. 18: Probability distribution of length of the spacer, and 
the length of the first bubble, at critical force g c = 0.82 and 
/3e = 3 (T = 1/3) for a disorder strength A = 5/3. A fitting 
function of the form f(m) = Ae~ am is used to fit the data. 
The best fit gives the value c = 0.72±0.01 for zipped segment 
and c = 0.31 ± 0.01 for the first bubble. 



with c b = 0.31 ± 0.01 for the bubble. 



VII. SUMMARY 

We obtained the exact phase boundary of unzipping 
of a polymer in 1 + 1 dimensions from a wall sepa- 
rating two regions of different affinity for the polymer. 
We developed a procedure to generate the phase dia- 
gram from an exact numerical procedure and for the 
pure case the numerical results do reproduce the exact 
results. The unzipping transition is characterized by a 
set of exponents describing the length dependence of the 
extensibility and the specific heat. Numerical results on 
the exponents also agree with the exact results, namely 



2,< 



l,d c = 1,. 



1. 



In presence of random impurities distributed indepen- 
dently, taken from uniform deviate of width A, a free 
polymer wanders more than in a pure medium. In ad- 
dition there is a free energy fluctuation that grows with 
the length. It is natural that there is a sample depen- 
dence of the unzipping process but the response under 
large force also has inner structure characteristic of the 
local pockets of favorable sites. However the steps one 
sees in the response averages out and one gets a smooth 
quenched averaged isotherm. The phase diagram is de- 
termined both for a hard wall and a soft wall. The ex- 
ponents d x = 1.5, </> = 0.5 are different from the pure. 
The behaviour of the specific heat and finer details of 
the sample dependence need further studies. Probabil- 
ity distributions of various relevant quantities are also 
obtained. 



APPENDIX A: PHASE BOUNDARY FOR PURE 
CASE 

In this appendix we give details to calculate the phase 
boundary for a pure problem. As stated earlier, we con- 
sider a more general case by assuming that there is a 
potential V(V > 0) at each lattice site on one side of 
the wall. The softwall and hardwall are then the limiting 
cases for V — > and V — > oo respectively of the problem. 
The recursion relation satisfied by the canonical partition 
function in the presence of potential is 



[Z N (x + 1) + Z N (x - l)]e-' 3y ,for x < 
Zn + i(x) = ^ Zn(x + 1) + Zn(x — 1), for x > 

[Z N (+1) +Z N (-l)]e 0e , foriE = 

(A-l) 

with the initial condition Zq(x) = e^ e S X) o. The generat- 
ing function for the partition function, 



G(z, x) = J2 z N Z N {x), (A-2) 

N 

can be taken to be of the form (ansatz) 

si i \ f for a; > . . „. 

G(z ' x) = \A'-^,forx<0 (A " 3) 

with A, A' and A to be determined, (z dependence of 
A, A' suppressed). The root test of convergence then tells 
us that the singularity closest to origin (on the positive 
real axis) in the complex z-plane determines the partition 
function for N — > oo. For finite N , a contour integration 
with deformation around the singularities of G(z,x) in 
the z-plane would yield Zm from Eqs. IA-2l and IA-31 
Using the ansatz Eq. IA-31 in Eq. IA-21 we get 

- = [A'-*" 1 + A'-^ 1 ] e-? v , for x < 0,(A-4a) 

— = [A* +1 + A*" 1 ], for x > 0, (A-4b) 

z 

- = [(A' + A)A+-1 e f)e , fon = 0, (A-4c) 

z z 

from which one obtains 



A = 
A' = 
and A = 



1 - VI - 4z 2 
2~z 



1 - Vl - 4z 2 e-^ v 
2ze-P v 
1 



1-(A' + A)ze^' 
The singularities coming from A and A' are 

1 A > 1 

2' 1 2cxp(-PV) 



(A-5a) 

(A-5b) 
(A-5c) 

(A-6) 
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A has the singularity 
1 



"-2 



1 



l + e-' 3V (l-2( 



(A-7) 



which depends on both the adsorption energy and the 
potential V. On taking the limit V — » 0, the problem re- 
duces to the softwall case. On the other hand the problem 
reduces to a hardwall case in the limit V — » oo. In these 
limits 22 becomes 



- (1 - e"^) 2 , for softwall 



(A- 



-y'l - (1 - 2e-' 3e ) 2 ,for hardwall 



The relevant partition function in fixed distance en- 
semble for large N is approximated by 



Z N (x) 



for x > 0, 



(A-9) 



via the contour integration method, with z 2 as the 
closest-to-origin singularity. The free energy is then 



(3F(x) = N\nz 2 -xlnX(z 2 ), 



(A-10) 



upto x-independent additive constant. The force re- 
quired to maintain the distance x is given by g = 
2dF/dx. A factor of 2 is needed as per our definition 
of unit length as the diagonal of a unit square of Fig. 
The phase boundary is then given by 



5c (T) = -2TlnA(z 2 ), 



g c (T)=T In 



1 + e-f> v (l - 2e~ 0e ) 
e -/3e(! _ e -pv e -p e ) 



(A-lla) 



(A-llb) 



which has been quoted in Eq. (|14a(l . On taking the 
limit V — > and V — * oo one gets respectively the phase 
boundary for softwall and hardwall. 



9c{T) = 



Tin [2e 0e - l],for softwall 
Tlnfe^-l], for hardwall 



(A-12) 



The zero force melting takes place at T c = oo for softwall 
and T c = e/ In 2 for the hardwall case. There is a nonzero 
T c for any V < oo. 

In the fixed force ensemble, apart from above men- 
tioned singularities, an additional force dependent singu- 
larity, Z3(/3g) = [2 cosh(/?<j)] _1 , comes from the generat- 
ing function (grand partition function) 



The phase boundary in the fixed force ensemble comes 
from equating two singularities z 2 — z$. It is 



g c {T) = 2Tcosh" 



1-1 



2it(l — vu) 
l + v(l- 2u) 



with 



e /3e , and v 



-0V 



(A- 14a) 
(A-14b) 



On taking the appropriate limits and with a little simpli- 
fication, one gets 



9o{T) 



2TtantT 1 [l - e'^] , for softwall 

2Ttanh _1 [l - 2e _/3e ],for hardwall 

(A- 14c) 

With some algebra one can check that Eqs. IA-14cl and 
IA-12l arc identical. 

On pulling the polymer on wrong side (x < i.e 
against the potential V), the phase diagram modifies 
as follows: In the fixed distance ensemble, the relevant 
free energy comes from the partition function Z^(x) ~ 
\'- x (z 2 )/z? +1 which gives the phase boundary 



g c {T) = 2T\n\'(z 2 ) 



(A-15) 



In the fixed force ensemble the force dependent sin- 
gularity, z 3 ((3g), of the generating function (Eq. IA-13|) 
modifies to z' 3 ((3g) = [2v cosh(/3g)] _1 which on equating 
with temperature dependent singularity, z 2 ((3e) gives 



9c(T) 



-2Tcosh" 



2it(l — vu) 
l + v(l~ 2u) 



(A-16) ' 

The low temperature behaviour of the phase boundary 
can be understood by a simple analysis. Since the wall 
is at x = 0, only the even monomers can be on the wall. 
For the soft wall case, this means the ground state has 
a degeneracy of 2 N for 2A^ steps because both sides of 
the wall are equally allowed. For small enough T, when 
a length 2m is stretched out by the force, the net change 
in the free energy would involve the loss of the bound 
state entropy and the gain in the energy due to stretching 
[Hi EH. This give F(m) = me - gm + Tmln2, from 
which the low temperature phase boundary comes out as 
9c (T) = e + T]n2. 

For any V > 0, the degeneracy is completely lifted as 
it is energetically not favourable to be on the negative 
x side of the wall. Consequently, there is no extra loss 
of entropy of the bound state at low temperatures. This 
explains the difference in the low temperature behavior of 
the phase boundary betweeen the softwall and the hard 
wall case (including V > 0). 



1. Finite size behavior 



N=0 



(A-13) Let us use the ansatz of Eq. IA-31 in Eq. I A- 1 31 for 

finite N, the sum over x goes from — N to N. Therefore 
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we shall have terms of the type [A(22)e^ s ] where Eq. 
IA-9l has been used. Noting that the critical force is given 
by Eq. lA-llal the N dependent term (that vanishes for 
N — > oo) is e N P( 9 ~ 9c \ One can then identify the finite 
size scaling variable as N(g — g c ) as quoted in Eq. HTH 



When the temperature exceeds T c (g) the polymer gets 
desorbed. In this case the thermodynamics is governed 
by the force dependent free energy iV _1 ^ r (T, g) = Tin 2:3 
where Z3 = l/[2cosh(g/T)] is the force dependent singu- 
larity. Therefore the fixed-force specific heat is 



APPENDIX B: SPECIFIC HEAT FOR THE PURE 
PROBLEM 



For a given force g, the temperature, T c (g), at which 
the transition from the adsorbed to the desorbed phase 
takes place can be calculated by using Eq (|A-12() . Below 
T c (g) the polymer is adsorbed on the wall and the ther- 
modynamic properties come from the free energy con- 
taining the binding term. The relevant partition function 
comes from z^ of Eq. IA-81 For the hard wall case, the 
corresponding free energy is N~ x T(T,g) = Th\Z2(T). 
The specific heat, which is the second derivative of the 
free energy is 



cm = t 



d 2 HT,g) 

dT 2 



e 2 exp(e/T) 



2T 2 (exp(e/T) - 1)' 



(B-E 



C g (T) = T 



d 2 T{T,g) 5 2 sech 2 ( 5 /T) 



dT 2 



T 2 



(B-2) 



There is a discontinuity at T c (g) but superposed on that 
there is a delta function for g ^ 0. A delta function at 
T c (g) shows that the transition is first order for g ^ 0. In 
the absence of force (g = 0) , only the discontinuity sur- 
vives and it is then a classical second order phase transi- 
tion. 
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